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Our understanding of localization in the integer quantum Hall effect is informed by a combination 
of semi-classical models and percolation theory. Motivated by the effect of correlations on classical 
percolation we study numerically electron localization in the lowest Landau level in the presence 
of a power-law correlated disorder potential. Careful comparisons between classical and quantum 
dynamics suggest that the extended Harris criterion is applicable in the quantum case. This leads to 
a prediction of new localization quantum critical points in integer quantum Hall systems with power- 
law correlated disorder potentials. We demonstrate the stability of these critical points to addition 
of competing short-range disorder potentials, and discuss possible experimental realizations. 



I. INTRODUCTION 

One of the hallmarks of the quantum Hall effect is the 
appearance of plateaus in the Hall resistance of a two- 
dimensional electron e function of applied mag- 
netic field. Transitions between plateaus have been asso- 
ciated with a quantum critical poin^. This critical point 
corresponds to a localization-delocalization transition for 
the electron's wave-function. Unlike the Anderson metal- 
insulator transition here the transition is between two in- 
sulating states: an extended wave function emerges at a 
single value of the magnetic field corresponding to the 
center of a Landau band. Experiments and numerical 
simulations have revealed many important properties of 
this critical poin^, such as the value of the correlation 
length exponent, Vq ~ 2.3. Attempts at formulating an 
analytically tractable theory of the transition have by 
and large failed. 

A very fruitful approach to unraveling the physics of 
plateau transitions has been the semi-classical descrip- 
tion of electron dynamics in the lowest Landau level^. 
Classically the electron's motion can be described as a 
fast cyclotron rotation in the plane accompanied by a 
slow E X B drift of the center of the circular orbit along 
lines of constant random potential, V(ii) =:const. Quan- 
tum mechanics adds tunneling and interference to the 
mix. 

The purpose of this work is to understand the rela- 
tion between classical and quantum dynamics of elec- 
trons in the lowest Landau level, in the presence of a 
random potential. In particular we investigate the effect 
of the fractal geometry of the classical orbits on the criti- 
cal properties of the quantum localization-delocalization 
transition, such as the value of Vq. One of our main re- 
sults is that we have identified a family of random poten- 
tials for which the classical electron orbits have continu- 
ously varying fractal properties, while the quantum crit- 



ical point remains unchanged. These random potentials 
are characterized by decaying power-law correlations in 
space. When the decay of correlations is not too fast we 
find new quantum critical points which are the quantum 
analogs of correlated percolation^. 

Previously, a real-space renormalization group treat- 
ment of correlated quantum percolation appeared in 
Ref. 1^. Preliminary results of our numerical approach 
were reported in Ref. |^ Here we provide a detailed ac- 
count of the numerical methods used to study the effect 
of power-law correlated disorder potentials on the inte- 
ger quantum Hall transition (IQHT), as well as theoret- 
ical arguments supporting them. We also describe new 
results on the stability of these quantum critical points 
to the presence of short-range correlated disorder. The 
stability analysis provides connections with experimentsii 
where these new quantum critical points might be ob- 
served. 

The work presented here addresses two different as- 
pects of the problem. On one side the motivation is pro- 
vided by understanding the role of disorder in determin- 
ing the transport properties of quantum Hall systems. 
On the other side, we are motivated by the possibility of 
extending correlated percolation results in the classical 
regime to the quantum regime. 



A. Disorder and the quantum Hall effect 

The role of disorder in the quantum Hall effect was 
recognized shortly after its discovery. From scaling the- 
ory it is known that electron wave-functions are localized 
by potential disorder in two-dimensions^. This remains 
true in the presence of a strong perpendicular magnetic 
field for states in the tails of the Landau bands, and leads 
to the observed plateaus in the Hall resistance^. Toward 
the center of a Landau band the localization length in- 
creases signaling the appearance of a delocalized state. 
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This results in the observed transition between plateaus 
in the Hall resistance and the peak in the longitudinal 
resistance. 

The precise nature of the disordered potential expe- 
rienced by electrons in quantum Hall systems remains 
poorly understood. The canonical picture (say, for GaAs 
heterostructures) is that ionized donors located between 
the two-dimensional electron gas and the sample surface 
are responsible for the random potential. Important ex- 
perimental progress has come from scanning techniques 
which provide a map of the electrostatic potential seen 
by electrons in a semiconductor heterostructurei^i. In- 
terestingly enough the experiment measures random po- 
tential fluctuations that extend over distances that are 
more than an order of magnitude larger than the typical 
spacing between the impurities and the electron gas, sug- 
gesting that the potential might be longer ranged that 
previously thought. Scanning techniques and the abil- 
ity to manipulate single atoms on surfaces also opens up 
the possibility of engineering the random potential so as 
to produce novel transport effects. For example, an en- 
hancement of the conductivity (compared to the Drude 
value) of a two-dimensional electron gas was observed in 
Ref. |3 due to the presence a random magnetic field. In 
this case the random field was produced by a thin film 
of rough magnetic material brought in close proximity to 
the electron gas. The fluctuations of the field were inher- 
ited from the height fluctuations of the films surface. 

Inspired by these experimental advances here we de- 
scribe how power-law correlated disorder can lead to 
novel critical behavior of the electron gas in the inte- 
ger quantum Hall setting. We also analyze the effect 
of competing disorders, which are necessarily present in 
semiconductor heterostructures, on the stability of these 
critical points. 



B. Quantum percolation 

The classical percolation problem provides one of the 
most intuitive examples of a critical point. Lattice sites 
are occupied with some probability p. Nearest neighbor- 
ing sites which are occupied form clusters with a typical 
size ^. As p is tuned to its critical value Pc, which depends 
on the lattice type, the correlation length ^ diverges as 
^ ^ \Pc ~ p\~'^'' ■ t'c is the correlation length exponent, 
and in two dimensions z/c = 4/3 independent of the lat- 
tice typaii. 

An analogous picture has been put forward to describe 
the localization transition in the integer quantum Hall 
system^. For a smooth random potential (one that varies 
little on the scale of the magnetic length) the electron 
wave- function is localized along the level lines. For ran- 
dom potential symmetric around = level lines away 
from the zero level are closed and have a typical size ^ 
which is identified with the localization length. As the 
electron's energy is tuned to the center of the Landau 
band the corresponding level line approaches V = and 



^ diverges. If the random potential has short range corre- 
lations in space this is precisely the classical percolation 
problemi^ and one would predict a localization length 
exponent z^^ = = 4/3. 

Experiments that measure a localization length expo- 
nent in the integer quantum Hall system find Vq ~ 2.3, 
signaling a break-down of the simple percolation picture. 
The essential physics that has been left out is one of elec- 
tron tunneling, which readily occurs at the saddle points 
of the random potential, as well as quantum interfer- 
ence. This was beautifully demonstrated by Chalker and 
Coddington who proposed a lattice model which takes 
into account these purely quantum effects, by describing 
electron dynamics as hopping from one saddle point to 
the next with scattering matrices associated with each 
saddloi^. The saddle points themselves occupy the ver- 
texes of a square lattice. Computer measurements of the 
localization exponent for this model yield values in agree- 
ment with the experimental results. Moreover, taking a 
classical limit of the network model, in which the scat- 
tering matrices become classical probabilities for the two 
possible outcomes of a scattering event at a vertex, leads 
to classical percolation and a localization length exponent 
of A/3—. This provides an intriguing connection between 
the critical points of classical and quantum percolation. 
Here we investigate the nature of this connection when 
power-law correlations are introduced in the random po- 
tential. In the classical case it was shown that this can 
lead to new critical points if the power-law decay is not 
too fast. We show that a similar effect occurs for the 
quantum version of percolation provided by electron lo- 
calization in the lowest Landau level. 

The paper is organized as follows. In section^]we dis- 
cuss the physics of classical and quantum electron mo- 
tion in a disordered potential in the presence of a strong 
magnetic field. We present scaling arguments that show 
how the value of the critical exponent v can be extracted 
from the time dependent mean square displacement (clas- 
sical motion) and the time dependent wave-packet spread 
(quantum motion). In section 11111 we describe in de- 
tail the numerical methods used to obtain the value of 
v in the classical and quantum setting. As evidence 
of the suitability of these methods, we present results 
for short-range correlated disorder potentials that are in 
good quantitative agreement with previous experimental 
and theoretical worksSiiiii. Section Hvl contains the nu- 
merical results obtained for classical and quantum elec- 
tron motion in power-law correlated disorder potentials. 
In section we propose a framework based on the clas- 
sical Harris criteriorki^ and its extensioni^ that provides 
the theoretical support for our numerical results. 

Finally, in section IVII we present results on classical 
and quantum localization in the presence of compet- 
ing short-range and power-law correlated disorders. In 
experiments^, which we believe have a good chance of 
observing the newly discovered quantum critical points, 
both types of disordered potentials arc present. The key 
question here is whether correlated quantum percolation 
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is stable to the addition of short-range correlated disor- 
der. 



II. CLASSICAL AND QUANTUM MOTION IN 
THE LOWEST LANDAU LEVEL 

The methods and tools used to describe the 
localization-delocalization transition in the lowest Lan- 
dau level (LLL) in the classical regime are considerably 
different from the ones used in the quantum regime. In 
this section, we review both approaches and present the 
arguments leading to the connection between the local- 
ization exponent v and electron dynamics. 

A. Classical motion 

The two-dimensional semi-classical motion of an elec- 
tron in a disorder potential under a strong constant mag- 
netic field (along the z-axis) , is described by the drift mo- 
tion of the electron's guiding center along equipotential 
lines: 

§ = ^Vnr)x£; (1) 

where r is the position vector for the guiding center in 
the (a;, y) plane. This equation is derived in the adiabatic 
approximation, i.e., under the assumption that the cy- 
clotron radius Ic = w/wc is much smaller than the typical 
distance over which the potential changesii. Therefore 
the localization-delocalization transition in the LLL in 
the classical regime, implies the study and characteriza- 
tion of the trajectories determined by this equation. 

Eq. was extensively studied in two different albeit 
related contexts. F. Eversi^used this equation in numer- 
ical studies of classical motion along the hulls of perco- 
lation clusters. He demonstrated interesting scaling be- 
havior of the time-dependent density-density correlation 
function close to the percolation critical point. 

V. Guraric tt alm^ arrived at Eq. from the classical 
limit of the Liouvillian equation of motion for electrons 
in the lowest Landau level of the integer quantum Hall 
regime. In both cases, the authors considered smooth 
random potentials with short-range correlations in space. 

One goal of this work is to study numerically properties 
of closed trajectories determined by Eq. where V{r) 
is power-law correlated in space. In particular, we focus 
on the scaling law, and its associated critical exponent v, 
which describes how the size of a closed trajectory grows 
as the energy approaches the critical value determined 
by the V{r) = equipotential. 

To make explicit the connection between v and the dy- 
namics of a particle moving in a two-dimensional random 
potential, let us briefly review the argument presented in 
Ref. For a fixed constant value of V{r) = Vb, there 
is a set of equipotentials associated with closed electron 



trajectories. These trajectories are the outer boundaries 
(hulls) of percolation clusters, with the occupation prob- 
ability p of the percolation problem determined by Vq- In 
particular, the critical value of Pc = 1/2 corresponds to 
Vb = in the units chosen. Several studie9i2ii242& have 
shown that the mean squared displacement between two 
points on a given trajectory Ar^ (r is the position vector 
for the particle on the trajectory), when averaged over 
all trajectories for fixed Vb, follows a diffusive law: 

{Ar'{t))v„ ^Dt (t « t*). (2) 

Here t* is a characteristic time that depends on Vq. 

Away from the percolation critical point particle tra- 
jectories are closed, and in the long-time limit the mean 
squared displacement reaches a constant value 

{Ar^{t))v„ = e [t » i*), (3) 

where ^(Vb) is the localization length. 

These considerations lead to a scaling form (at fixed 
Vb) for the mean squared displacement 

(Ar^(t))yo = ^^/(^) ' (4) 

where / is a scaling function. The classical result from 
percolation theorySS for the correlation length, ^ oc |p — 
Pc\~'^'' ^ leads to the expression ^ oc V^^", which in tm'n 
implies 

{Av\t))y, = Dtfiy^}i. (5) 

Finally, by averaging this equation over all values of 
Vb one obtains the scaling relation 



(Ar2(t)) ^ t', (6) 

where 9 = \ — \/2i'c is the anomalous diffusion exponent. 
From simulations (Ar2(t)} can be computed which leads 
to a value for 9 and, via the scaling relation just derived, 
to a value for the localization length exponent Vc in the 
classical regime. 

B. Quantum motion 

One approach to study the localization-delocalization 
transition in disordered quantum systems is based on 
the different contributions that extended and localized 
electron wave- functions make to the frequency depen- 
dent electrical conductivitySi. Localized states appear in 
the electrical conductivity through the retarded density- 
density Green's function, a quantity that has an intuitive 
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interpretation when written in terms of wave- functions. 
When the electron is represented by a locahzed wave- 
packet at position r and time t ~ Q, the Green's function 
gives the value of the wave-packet spread at time t. Since 
this wave-packet describes the probability of finding the 
electron at a given position and time, its spread is a mea- 
sure of the uncertainty on the electron's position, and as 
such, an indicator of localization. 

In most studies of localization, the Green's function 
is calculated in the momentum-frequency domain and 
studied in detail in the limit of large momenta/small 
frequenc}i2*2S. The purpose of this section is to show 
explicitly how the critical exponent can be extracted 
from the Green's function using instead the real time do- 
main, as proposed by Sinova, Meden and Girvin^i. 

Within this approach, one studies the disorder aver- 
aged density-density correlation function projected onto 
the LLL, i.e., the focus is on the unconstrained spec- 
tral function. As pointed out in Reffioi. the applicability 
of the method rests upon the assumption that delocal- 
ized states in the IQH transition are isolated from each 
other and located at the centers of the Landau bands. 
Thus, when the density-density correlation function is re- 
stricted to one Landau level (the lowest being considered 
for simplicity), the contribution of the (single) delocal- 
ized state in that level dominates the sum over states 
in the small (<?,w) limit. Equivalently, the integral of 
the spectral function over all energies is dominated by 
the contribution from the critical energy in that limit. 
For our purposes, the main advantage of this approach is 
that it clearly spells out a numerical scheme for studying 
the effect of power-law correlated disorder potentials on 
localization in the lowest Landau level. 

In order to make connections with electron motion in 
the classical regime and also to fix the notation used in 
the rest of the paper, we review some of the main points 
of the approach. Consider 2d spinlcss electrons (the elec- 
tron spin is fixed by the magnetic field and we assume 
that spin-flipping interactions are not allowed) in the x- 
y plane, under the combined effects of a magnetic field 
B = Bz, perpendicular to the plane and a random po- 
tential V{r) due to the presence of impurities. 

The electron-impurity interaction is given by: 



with the projected one-particle density operator taking 
the form: 



(9) 



with C, the position vector for the guiding center of a 
one-particle orbit, defined as24: 



eB " 



(10) 



Here n = p + |A and (x, y) are the canonical momen- 
tum and position operators for one particle. 

The formalism used to project the density operator 
Pq = e^**! "" onto the LLL was developed in Refs. I25ll2^ 
There it was shown that the commutation relation for the 
projected density operator obeys a closed algebra (the 
magnetic translation operator algebra): 



[Pk , Pq] = 2i sin i^-^k xqj e ^'^^'^ pk-(-q, (11) 

with k X q = (k x q) • i. This property implies that the 
equation of motion for the density operator 



(12) 



is closed and can, in principle, be solved exactly. 

Following the same line of argument, one can show 
that similar results hold also for the correlation function 
of the projected density operator: 



G(r,t) = rr(p(r,t)p(O,0)) 



(13) 



note that the trace extends over states in the lowest Lan- 
dau level only. The density-density correlation function 
G in k-space satisfies the equationS^: 



h—G{k,t) = ^ 2 i sin ( ^fc X g] F(k- q 



kPk 



(7) 



exp[-|(fc2_k.q)]G(q,i) (14) 



where pk 



„?k-r 



is the one-particle density operator 



and Vk is the Fourier transform of the disorder poten- 
tial. At high enough magnetic fields (or low enough tem- 
peratures), the kinetic energy is quenched to the value 
of the lowest Landau band and the sum is reduced to a 
sum over states in the LLL. Let H and p be the Hamil- 
tonian and one-particle density operator projected onto 
the LLL. Then, 



H 



k Pk, 



(8) 



which follows from Eq. This is a Schroedingcr-like 

equation where the effective Hamiltonian is the Liouvil- 
lian matrix: 



/:kq = 2zsin(^|fcx9^ T/(k-q) exp[-|(fc2 _ k • q)] 

(15) 

(The name Liouvillian is used because of the analogy 
with the Liouvillian operator in classical mechanics.) 

An interesting feature of Eq. H14|l is that for a fixed 
value of the magnetic field, the magnetic length €c {^t = 
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vanishes in the classical limit h ^ 0. In units of 
^ = 1, Eq. (ini) givesia: 

^G(r,t)=e,,9,l/(r)a,G(r,t). (16) 
A solution of this equation is G(r, t) = 5{r—r{t)) where 



— ^VV{x,y)xz. (17) 

In other words, Eq. H16() describes the classical motion 
of the electron's guiding center along equipotentials of 
V{rp.. 

As argued above, the Liouvillian driving the quantum 
dynamics of an electron in the LLL, contains information 
about localization and the value of the critical exponent 
Vq. The following scaling argument provides a numer- 
ical algorithm for extracting the localization exponent 
from the time dependent projected density-density cor- 
relation function. (We would like to draw attention to 
a close analogy between the present argument and the 
ideas presented in Sec. Ill Ah 

Let us start with Eq. (jH)) , and assume that all eigenval- 
ues and eigenstates, localized and delocalized, are known. 
Assume next that the disorder potential breaks com- 
pletely the degeneracy of the LLL, and that there is only 
one extended state at energy E = 0. We construct a 
wavepacket at time t ~ using only localized states with 
energies in the A neighborhood E' ^ 0, and then we let 
the wave-packet evolve in time. The state-vector corre- 
sponding to such wave-packet is: 

E'+A 

where \i > are eigenvectors of the Hamiltonian projected 
onto the LLL (see Eq. ©). 

At very short times, the spread of the wave-packet as a 
function of time is ballistic (see Appendix for derivation) . 
However, after some crossover time and before the local- 
ization length corresponding to the energy E' is reached, 
the spread of the wave-packet averaged over disorder re- 
alizations is diffusive: 



< Ar^{t) >E' ^ Dt (t < t*{E')) . (19) 

At times much longer than t* ~ C^, the localization 
length is reached and 

<Ar^t) >E' = e{E') it » t*iE')) . (20) 

In close analogy with the classical argument, a scal- 
ing form for the average dispersion of the wave-packet 
follows: 



<ArHt)>E'^Dtr(^-^'^ , (21) 

where /' is the appropriate scaling function. Now if a 
localized wave-packet at time t = is constructed with 
all the eigenstates (including the lone delocalized state 
at the center of the LLL) , its wave- vector can be written 
as: 

m)>= E \^E'{t)> . (22) 

E'^-oo 

Using this decomposition the dispersion of \ip{t) > can 
be expressed as 

^■E'=o 

< Ar^t) > ^ / dE' < Ar2(i) >e- g{E') , (23) 

J E'^-oo 

where g{E') is the density of states; the contribution 
of the cross terms can be neglected assuming that the 
extended state is non-degenerate. Specifically, it can be 
shown that the cross terms are a correction to the leading 
behavior if there is only one extended state at the critical 
energy [E = 0). Replacing < Ar'^{t) >e' by Eq. 12}, 
using the fact that the density of states is non-singular 
and that ^ is related to E' by £,[E') ^ [E')-"" , we obtain: 

< Ar^t) > = /" dE' Dtg{E')f'{-^) (24) 
Finally, an appropriate rescaling gives: 



<Ar2(t)>cxr; e = l~l/2vq, (25) 

and we conclude, as in the classical case, that the spread 
of the wave-packet is sub-diffusive with an anomalous 
diffusion exponent 6. 

Therefore, by computing < Ar'^it) >, the value of v 
can be obtained. An alternative derivation of this re- 
sult has been proposed in Ref. |23, where the energy in- 
tegrated Liouvillian propagator is analyzed in the finite 
(q, w) regime and the limits of g ^ and small ui are 
taken. 



III. SHORT RANGE DISORDER POTENTIAL 

We introduce next the numerical procedures used to 
calculate the value of the critical exponent v. We do 
this first in the case of a short-range correlated disorder 
potential, where we test the effectiveness of these numer- 
ical techniques by re-deriving previously known results. 
We begin by introducing the method used for analyzing 
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classical localization, where the connection to percola- 
tion theory provides a testing ground for our numerics. 
In subsection IIII Bl we describe how the Liouvillian ap- 
proach is implemented following the method introduced 
in Ref. to analyze quantum localization. The numer- 
ical values thus obtained compare favorably to experi- 
mental results and previous numerical calculations. 



A. Classical motion 

To study the classical motion of the electron's guiding 
center along the equipotentials of ^(r) we make use of 
a lattice model, used previously to compute geometrical 
exponents for contour loops of rough interfacesSS. 

Consider a square lattice C with N sites and periodic 
boundary conditions. Each site in the lattice has assigned 
to it a random number that represents the value of the 
disorder potential V{y) at that point in space. Its dual 
lattice £* is the set of points that describe the positions 
of the electron in this potential landscape. Thus, an elec- 
tron originally located at a site Tq of C* will move along 
a path that joins points on this lattice. The trajecto- 
ries are contour loops of the disorder potential. They are 
generated by first randomly selecting the initial position 
ro from the set of points in C* . The choice of the ini- 
tial position of the trajectory also fixes the value of the 
level Vq. Vb is computed as the average of the values for 
V{y) at the four vertexes (sites in C) that surround the 
initial position Tq. Once Vo is determined, all the lat- 
tice sites of C can be labeled by + or — , depending on 
whether they are at a potential higher or lower than Vq. 
The numerical algorithm updates the electron's position 
by selecting the bond on the dual lattice that connects 
the point Tq with one of its nearest neighbor points, with 
the additional property that it crosses a H — bond on the 
lattice C. With this construction, the contour loop is a 
directed walk along the bonds of the dual lattice £* , that 
separates potentials lying above (below) Vq on the inside 
(outside) of the loop (see Fig.QJ. 

At each step, the square of the displacement vector 
between the current and original positions on the dual 
lattice is calculated. The procedure is repeated for dif- 
ferent starting points rg, and different disorder realiza- 
tions. All the calculated quantities are averaged over all 
trajectories and over all disorder realizations. 

A technical subtlety of this procedure is the existence 

of saddle-point plaquettes, where there are H 1 — signs 

cyclically around them. In these cases, two contours (four 
links) meet at the center and it is necessary to add an 
additional rule to resolve the connectivity, so as to con- 
vert this pattern into two 90° turns that are not quite 
touching. A physically sensible rule makes use of the av- 
erage potential Vpiaq of the four potentials around the 
saddle-point plaquette. If Vpiaq < Vq, the center of the 
plaquette is at a lower potential than the contour loop 
and the connectivity is resolved by having the -I- sites in- 
side the 90° turns. In the opposite case, Vpiaq > Vq, the 
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FIG. 1: Classical trajectory of the electron's guiding cen- 
ter, in the presence of a random potential whose values are 
indicated. Vb is the average potential around the plaquette 
surrounding the initial point of the trajectory. The resolu- 
tion of a saddle point is determined by the average potential 
around the plaquette surrounding the saddle. 



-I- sites are outside of the 90° turns (see Fig.P). 

As in Ref. we find from simulations that the average 
time for an electron to traverse a certain distance on C* 
is proportional to the distance. This observation justifies 
using the number of steps along C* to measure the time 
elapsed. In this way time is rendered dimensionless. 

The computations were carried out on two-dimensional 
square lattices with system sizes N = 256, 512 and 1024. 
The values for the short-range correlated potential were 
chosen independently from site to site from a uniform 
distribution over the range [—0.5,0.5]. This choice leads 
toV — for the value of the critical equipotential energy. 
Fig. m shows results for a 1024 x 1024 lattice obtained 
by averaging over 3 x 10^ trajectories and 1000 different 
disorder realizations. 



As discussed above, the diffusive (slope 9 = 1), short- 
time {t < t* ^ 10), behavior is followed by a crossover 
to a critical regime where the value of the slope 6* is a 
measure of the critical exponent Vc', 0=1 — l/2vc- By 
fitting the critical region to a line on a log-log graph we 
find 



1.25 ±0.05 



(26) 



where the error bars rcfiect the uncertainty associated 
with choosing the critical region. 

This value compares well with the exact result from 
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FIG. 2: Averaged mean square displacement for a classical 
particle drifting along equipotentials of a random short-range 
correlated potential. Distance is measured in units of lattice 
spacing and time in number of lattice steps. System size is 
1024 X 1024. V is related to the slope in the critical region Q, 
by e = 1 - l/2i^. 



The electron density operator p{ki,k2) = 
exp(2'Ki(kix + k2y)) projected onto the LLL us- 
ing this basis, is a matrix of the form: 

p(fci,fc2) = cxp (^-M_LM^^ i(fci,fc2) (28) 

where ki, k2 are integers ((fci, fc2) (0, 0)), and the ma- 
trix L(ki , fc2) is given in terms of two unitary unimodular 
NxN matrices: 

L(fcl,fc2) = £'=1^2/2 jfci ^fe2 

/ 1 0\ 



/ = dzag(l,e,...,e^-i) 

with e = exp(2^); note that ft. is a cyclic permutation 
matrix. These matrices satisfy: h f = e f h and = 
jN _ ^ rjij-^g explicit expression for the matrix elements 
of L(fci, fc2) is: 



percolation theory = 4/3—. The critical regime 
extends until finite size effects become dominant and 
< Ar^ > saturates, as observed in Fig.[21for t>2x 10^. 



B. Quantum motion 



modN 



(29) 



Because the basis of states is restricted to the LLL, the 
kinetic energy is a constant and the projected Hamilto- 
nian, reduces to: 



Following the ideas presented in Sec lII Bl we studied 
numerically the role played by quantum effects on elec- 
tron dynamics in the lowest Landau level and in the pres- 
ence of a short-range correlated random pot ential. 

We use the approach proposed in Ref. [l^ and use the 
eigenstates of the Hamiltonian for an electron on a torus 
geometry as the basis of states for the LLL. The cor- 
responding wave-functions written in the Landau gauge 
are: 



cxp (27r(a; + iy){N in + a)- 

_m— — oo 

{Nm + af 



N 



(27) 



where m takes integer values, is the number of flux 
quanta through the torus, and the index a goes from 
to — 1, labeling the A^ states in the LLL. These wave- 
functions are periodic functions in the interval (x, y) G 
[0,1) X [0,1), and are centered along narrow strips (of 
width ^c) around the lines x = a/N] an example is shown 
in Fig.^Ji). In Eq. H27|) (x, y) are dimcnsionlcss variables 
expressed in units in which the magnetic length is tc= 
l/27rAf. 



H = J2 y{-ki,-k2)p{ki,k2) 



(30) 



^ Vi-k,,-k2)cxp{^'^^^^7T] L{kuk2, 



fci ,k2 



2N 



where the sum runs over all integer values of (fci, /s2), and 
1^(0, 0) = (this choice amounts to fixing the critical 
energy at zero). 

In this basis, the Hamiltonian is a random NxN ma- 
trix that can be diagonalized exactly. Notice that the 
projection operation introduces the exponential factor in 
Eq. that effectively reduces the amplitude of the 

large (fci, k2) Fourier components of the potential (it acts 
as a soft large- momentum cutoff). The presence of this 
factor seems to indicate that the transition is driven by 
the small momenta components of the disorder poten- 
tial. It is important to remark, that even when a nu- 
merical calculation procedure fails (due to large but fi- 
nite precision) to detect the contribution coming from 
momenta components higher than a maximum value for 
{kf + fc|)/2A^ (numerically, the Hamiltonian matrix be- 
comes a banded-matrix), it is incorrect to replace the 
exponential factor by a hard-cutoff, since the commu- 
tation relations of the projected density operator— are 
violated. (In numerical calculations, this translate into 
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having results that are extremely sensitive to the value 
of the hard-cutofF chosen.) 

The values for the Fourier components of the disorder 
potential in Eq. are chosen randomly from a flat 
distribution over the interval [—0.5,0.5] in units of HlUc- 

We diagonalize the Hamiltonian to obtain all N eigen- 
values and eigenfunctions. Figure 13 shows examples of 
a localized state at the tail of the Landau band, and a 
delocalized one close to the band center, obtained with a 
basis of = 1000 states. 




FIG. 3: Squared amplitude of a localized eigenstate with 
energy E = 9.6 (top), and a delocalized one with E = 0.011 
(bottom), in the lowest Landau level with degeneracy A'' — 
1000. The disorder potential is short ranged. 

A localized wave-packet is constructed with all the 
eigenstates for a given disorder realization and, as it 
evolves, its spread along the x-direction is computed as 
a function of time. The wave-packets are chosen from 
among the basis states, Eq. (|27|l . which are localized 
in the x-direction and completely spread out in the y- 
direction (see Fig. 0Jl). Thus, the contribution to the 
total spread from the y-direction is a constant propor- 
tional to the system size. Fig0]shows a set of 'snapshots' 
for the evolution of the wave-packet initially localized at 
X = 0.5. 

The procedure is repeated for different initial positions 



for the localized wave-packet and different disorder real- 
izations. The spread is averaged over all initial positions 
and disorder realizations. A typical result for a basis of 
1000 states and some 1000 disorder realizations is plotted 
in Fig|5| 

As in the classical case, three regimes can be identified. 
For long enough times (t > 10 in FiglS)) the spread of 
the wave-packet reaches a constant value indicating the 
influence of the finite system size. The critical region 
corresponds to intermediate times (0.4 <t< 10), and the 
value of Vq ("q" for quantum) can be related to the slope 
of the line, which is the anomalous diffusion exponent 9 
(see Eq. (|25|l ). As seen in the inset to Fig|5l the slope 
in the critical region shows a strong dependence on the 
system size. To take into account these finite size effects, 
we determined the slope for systems sizes ranging from 
A^ ~ 200 to A^ = 1500 and compute the value for an 
infinite sized system by linear extrapolation to obtain 
i^q = 2.33 ±0.09. 

Our result compares favorably to the values measured 
in experiments {ly = 2.3 ± 0.1)^ and previous numerical 
simulations (ly ~ 2.35 ± 0.05)^. This provides an impor- 
tant check on the numerical method. 

Our data shows that at very short times the slope 9^2 
which implies ballistic spreading of the wave-packet. In 
the Appendix we present a calculation based on pertur- 
bation theory that shows that this is a result of averag- 
ing the evolution equation over random disorder, at short 
enough times. 

IV. POWER-LAW CORRELATED DISORDER 
POTENTIAL 

In this section we address the central issue raised by 
the present work: how a change in the properties of the 
probability distribution for the disorder potential affect 
the nature of localization in the lowest Landau level. The 
numerical procedure outlined in previous sections allows 
us to investigate this point in a straightforward manner 
for both the classical and the quantum regime. Basically, 
it amounts to an appropriate modification of the proba- 
bility distribution for the Fourier components for the ran- 
dom potential in Eq. H31() . Since we are interested in the 
effect of power-law correlations on localization properties 
of the Hamiltonian in Eq. H31|l . the Fourier components 
are chosen independently with a Gaussian distribution 
with zero mean, and variance 

|l/(k)|2 oc l/|k|2-" . (31) 

Inverse Fourier transforming to real-space leads to power- 
law correlations 



y(r)F(0) - l/|r|" . (32) 

The constant in Eq. H31|l is fixed so that the vari- 
ance of V{r) is normalized to one. Long-range correla- 
tions strongly modify the real-space configuration of the 
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FIG. 4: Evolution of a wave-packet made of all eigenstates that is initially localized nt x — 0.5. The squared amplitude of the 
wave-packet is plotted as a function of the real space coordinates (x, y) at four different times. 



equipotential lines as can be seen in Fig. |^ where a com- 
parison between a short-range and a power-law correlated 
potential with a = 0.5 is given. 



A. Classical motion 

In the classical limit of the IQH transition, the effect 
of power-law correlations in the disorder potential can 
be described in terms of a purely geometric effect: they 
change the fractal geometry of the equipotential contours 
modifying the electron's path, and correspondingly, the 
anomalous diffusion law for its mean squared displace- 
ment. 

To study numerically this regime, we repeated the pro- 
cedure described in Sec lIII Al and computed {Ar^{t)) for 
values of a ranging from 0.33 to 1.8. Fig. [7| shows the 
curves obtained for a system of 512 x 512 lattice sites. 
As in the case of short-range correlations, three regimes 
can be identified with the onset of each of them depend- 
ing on the value of a. From the figure, a qualitative 



change in the spread of the electrons position in time 
is observed. Namely, in the critical regime (typically for 
30 < t < 1100), the value of the slope becomes dependent 
on a when a < 1.5. As the value of a decreases (correla- 
tions increase) below a critical value of a* ~ 1.5 ('c' is for 
'classical'), the slope 9 in the critical region increases and 
h' becomes an increasing function of a tending to infinity 
as a approaches zero. 



B. Quantum motion 

In contrast to the classical regime, the effect of power- 
law correlations in the disorder potential on the prop- 
erties of quantum localization, does not have a simple 
geometrical interpretation. In particular, we find that 
changing the fractal properties of the equipotentials does 
not always lead to a change in the quantum critical ex- 
ponents. 

The numerical procedure introduced in Sec. IIII Bl is 
easily extended to the present case allowing us to study 
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FIG. 5; Spread of wave-packet as a function of time, in the 
presence of short-range correlated disorder potential with a 
basis of A*' = 1000 states. The spread at intermediate times 
is sub-diffusive. The inset shows the anomalous diffusion ex- 
ponent 6^ as a function of basis size A''. The intercept gives 
the infinite-system-size value of Q which leads to the quoted 
value of the localization exponent (i' = 1/(2 — 26)). 



the problem in detail. The Fourier components of the 
random potential are taken from the same distribution 
as in the classical case and the corresponding quantum 
Hamiltonian (Eq. (|31|l ) is diagonalized. Then, < Ax^ > 
is computed as a function of time for values of a ranging 
from 0.1 to 1.9. A typical set of results obtained with a 
basis of 1000 states is shown in Fig. together with the 
curve obtained for a short-range correlated potential. As 
in the short-range case, the exponent Q was calculated 
for different systems sizes ranging from TV = 300 to iV = 
1000. 

Like in the classical case, a qualitative difference in 
the slopes of the critical regime develops as a changes. 
A careful comparison among the different curves shows 
that the value of the slope in the critical region starts 
to increase below a critical value of a* ~ 0.85. A de- 
tailed quantitative analysis of these results together with 
a comparison to the classical case is the subject of the 
following section. 



V. EXTENDED HARRIS CRITERION 
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FIG. 6: Contour plot of a typical realization of a short-range 
correlated potential (top) and a power-law correlated one with 
exponent a = 0.5. 

quantity which couples to the disorder, within a volume 
set by the correlation length ^, do not grow faster than 
its average value, as the transition is approached and 
^ — > oo. It was proved rigorously and generalized to a 
wider class of systems if an appropriate definition for the 
correlation length is adopted^. An extension of the cri- 
terion to power-law correlated disorder was proposed by 
Weinrib and Halperirki&. 

In the case of two-dimensional percolation with power- 
law correlations in the occupation probabilities, with 
the power-law exponent a, the extended Harris criterion 
states: 



The effect of a random disorder potential on the prop- 
erties of continuous phase transitions in classical systems 
is summarized by the Harris criterioni^. Basically, the 
criterion states that critical exponents of the disordered 
and the corresponding clean system remain equal as long 
as the value for the correlation length exponent v satis- 
fies: di/ — 2 > 0, where d is the dimensionality of the clas- 
sical system under consideration. The criterion is derived 
by requiring that the fluctuations of the thermodynamic 



a > a* = 3/2 - 

a < v,^- . (33) 

a 

As a numerical check on Eq. (|33|l . we computed the 
values of the anomalous diffusion exponent 9 for classical 
electron motion and plotted them as a function of a in 
Fig. El The values of 9 were obtained from fitting the data 
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FIG. 7: Classical motion: averaged mean square displacement 
as a function of time for 512 x 512 system size. Distances are 
measured in units of lattice spacing and time in number of 
lattice steps. 
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FIG. 9: Classical regime: exponent S as a function of power- 
law correlated disorder exponent a. The horizontal line rep- 
resents the theoretical value of 8 — 0.625 corresponding to the 
value of i/c = 4/3 obtained from 'classical' percolation theory. 
The linear dependence of6 = l— Q/4is the prediction of the 
extended Harris criterion, Eq. 1331 
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FIG. 8: < Ax^ > as a function of time for various values 
of the power-law correlated disorder exponent. Anomalous 
diffusion appears at intermediate times. The top four curves 
have been multiplied by a factor of 10 for clarity. 



following the same procedure as in the classical case. The 
values plotted arc the extrapolation corresponding to an 
infinite system size. We see that the numerical results 
presented in Fig. IIUI are in good agreement with the ex- 
tended Harris criterion which now reads: 



a > ce., 



a < a„ 



0.85 
0.85 



2.33 
2 

a 



(34) 



Support for the validity of the extended Harris cri- 
terion for quantum critical points was also provided by 
studies of random quantum magneta^i. 

In the case of quantum percolation (i.e., localization 
in the lowest Landau level) the extended Harris criterion 
can be argued for in close analogy to the classical casein. 
Namely, consider the average size of the fluctuations of 
the electron's energy in an area set by the localization 
length ^. Like the electron's energy, it is also determined 
by the random disorder potential, and can be approxi- 
mated as: 



in Fig. 13 in the critical region, to a line using a least- 
squares method; 9 is the slope of the line. We see that 
the data is well described by the theoretical prediction: 
above a* = 3/2 the anomalous diffusion constant 9 = 
1 — l/2:^c = 5/8 is constant, while below this value it 
varies as 9 = 1 — a/4. Similar results were obtained 
previously by Prakash et al^. 

Analogous results hold true for the quantum case. In 
Fig.^Jwe plot the anomalous diffusion exponent 9, which 
describes the spreading of an initially localized wave- 
packet in the lowest Landau level, as a function of a. 
The exponent 9 was computed using the data in Fig. |S1 



^2 « (^1^ d^xV{^)^ 

^ - ( d^a-ixj-" (35) 

where in deriving the second equation we have made use 
of < V{x)V{x') >~ l/|x - x'l". 

For the quantum critical point associated with local- 
ization in the lowest Landau level to be stable to the 
introduction of power-law correlations in the random po- 
tential, the fluctuations in the electron's energy should be 
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FIG. 10: Quantum regime: anomalous diffusion exponent 
9 = 1 — l/2vq as a function of the exponent a, which deter- 
mines the power-law correlations of the disorder. The lines 
are predictions of the extended Harris criterion (Eq. 1341 ') us- 
ing the value Vq = 2.33 for the quantum localization length 
exponent, in the case of short-range correlated disorder. 



small when compared to its energy E. Using Eq. (|35|l . 
and the scaling relation E ^ ^^^/'^ we arrive at the esti- 
mates 





a 


> 2 




E^'^-^lniE)-" 


a 


= 2 


(36) 


j^av — 2 


a 


< 2 . 





E^ 
E^ 



The quantum percolation critical point is expected to 
be unaffected by the introduction of power-law correla- 
tions if E'^/E'^ when E 0, or, equivalently, when 
^ — > oo. Since ly > I this is always the case when a > 2. 

However for a < 2 there are two regimes to consider. 
When 2 > a > 2/iy, E'^/E^ and the localization 
critical point due to short-ranged potential disorder is 
again stable to introduction of power-law correlations in 
the random potential. In contrast, for values of a < 2/i^ 
a power-law correlated potential produces large fluctua- 
tions in the energy thus destabilizing the critical point. 

The usual expectation is that this relevant perturba- 
tion leads to a new quantum critical point characterized 
by a new value of v. This is confirmed by our numerical 
results. 

An important consequence of the present result regards 
the validity of the expression put forward in Ref . 0, that 
relates the values of Vq and as follows: 



1 



(37) 



According to the analysis presented above, a quan- 
tum system with a long-range correlated potential has 
the same critical exponent as one with short-range corre- 
lations as long as a > 0.85. However, a classical system 
with potential correlations with a value for a in the re- 
gion 0.85 < a < 1.5, has a continuously varying critical 



exponent ly according to Eq. H33|) . For instance, a value 
of a = 1 gives i^q = 2.33 while i^c — 2.0. This analysis 
suggests that the relation given by Eq. (pITjl (whatever 
the argument used to support it) is inconsistent with the 
Harris criterion. 



VI. COMPETING DISORDERS 

In this section we undertake a numerical study of the 
effect of competing disorders on localization in the LLL. 
We investigate first the classical regime, which is in the 
universality class of correlated percolation. Analytical 
calculations^ have shown that in these systems the clas- 
sical critical point generated by a long-range correlated 
potential is stable against perturbations produced by a 
short-range correlated potential. As shown below, our 
numerical results in this regime, provide support for this 
picture. 

In the quantum case we investigate the effect of com- 
peting power-law correlated and short-range correlated 
disorders on localization. Numerical studies of this par- 
ticular case are relevant for experiments attempting to 
measure properties of the correlated quantum percola- 
tion critical point. Namely, since short-range correlated 
potentials are always present in experimental setups, it is 
important to understand quantitatively how their pres- 
ence affects the values of the new localization length ex- 
ponents described in the previous sections. 

To proceed with the appropriate numerical model we 
first establish the following conventions. The total dis- 
order potential is defined as a weighted sum of two po- 
tentials: one with short-range correlations ^^(r) and one 
with long-range correlations Vi{r). The units chosen are 
such that the mean values of short- and long-range corre- 
lated potentials are set equal to zero while their variances 
are normalized to one. In this way, the maximum am- 
plitudes of both potentials are comparable. In order to 
vary the relative strength of each potential we introduce 
a parameter k with values in the range < k < 1 such 
that the total potential is: 



Veffir) = KVs{r) + il-^)Vi{r) . 



(38) 



Notice that when k = the total potential has power- 
law correlations only, while the short-ranged part is in- 
creasingly dominant as the value of k increases toward 1 . 
This definition ensures that the variance of the total po- 
tential Veff remains fixed as k changes and (l/yy(r)} = 1. 

The numerics were performed with short-range corre- 
lated potentials constructed from random Fourier com- 
ponents I^(k) sampled from a imiform distribution, and 
from a Gaussian distribution with variance ^ |k|"^^ with 
a > a*. In real space this leads to delta-correlations 
and power-law correlations with exponent a, respectively. 
Long-range correlated potentials were generated with 
power-law correlations with exponent /3 < a*, for which 
the extended Harris criterion predicts a new localization 
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TABLE I: Statistical properties of the short-range and long- 
range correlated potentials used in section Wl\ 

Veffir) = /cK(r) + (l-/c) l/^r) 

{V4r)) = mir)) = 

{K'(r)> = iV^r)) = 1 

{Vs{r)Vs{r')) cx |r-r'|-", <5(r - r') 

(Vi{r)Vi{r')) cx |r-r'|-'' 



length exponent v = 2//3. The properties of the total 
potential and its components is summarized in Table 

A note on the evaluation of errors: as it will be shown 
in the following sections, the effect of tuning k 7^ de- 
pends on system size for both classical and quantum 
regimes. Hence, it was not possible to carry out the finite- 
size study similar to the one done previously, where there 
is only one type of disorder present. As a consequence, 
the errors have been estimated from the largest system- 
atic error involved in the procedure, which appears in the 
determination of the slope of the critical region. 



the value of 9 (and hence of ly) remains closer to the 
pure long-range value up to fairly large values of k 
(roughly 55% — 60%). The increasing sharpness of the 
crossover with increasing system sizes suggests a very 
sharp crossover in the thermodynamic limit. In the lan- 
guage of rcnormalization group, this is an indication of 
the stability of the long-range critical point when per- 
turbed with a short-range correlated potential. Notice, 
as it was pointed out above, that this crossover is strongly 
dependent on the system size, rendering useless the finite- 
size scaling analysis used previously. An interesting effect 
is observed in the raw data that points to the way the 
crossovers occur in finite size systems. Between the initial 
crossover from diffusive to sub-diffusive motion and the 
final crossover from sub-diffusive motion to saturation, 
there arc two clearly distinguishable regimes. At the ear- 
lier times, the slope is determined by the short-range 
correlated potential. Only at latest times (and before 
saturation effects dominate), a new slope appears that 
corresponds to the value of 9 determined by the expo- 
nent of the power-law (long-range) correlated potential. 



A. Classical Regime 

Computations were carried out for system sizes N = 
256, 512 and N = 1024. The values of the anomalous dif- 
fusion exponent 9 were extracted from the critical region 
of the spread of the classical position of the electron, as 
in section IIV Al Values of 9 were obtained for a range 
of values of the parameter k. Fig. ^2 shows result that 
correspond to a combination of a power-law correlated 
potential with exponent (3 = 0.5 (for the long-range part) 
and a = 1.8 (for the short-range piece). 




FIG. 11: Classical regime: slope as a function of param- 
eter K. K = 0(1) indicates a pure long-range (short-range) 
correlated potential. 

The figure shows that, as the system size increases 



B. Quantum Regime 

In contrast to the classical regime, up to date, there 
is no experimental or theoretical work to our knowledge 
that attempted to investigate the influence of mixed long- 
range and short-range correlated disorder potentials on 
the integer quantum Hall transitions. The numerical 
strategy used in section FlV Bl provides us with a frame- 
work to study this situation and allows us to examine the 
stability of newly found quantum critical points. 

Computations for quantum systems were carried out 
with a basis of iV = 500, 700 and N = 1000 states. The 
values of the slope 9 were extracted from the critical re- 
gion and plotted as a function of the parameter k, as in 
the classical regime. The results shown in Fig. ^] were 
obtained with power-law correlated potentials Vs and Vi 
with exponents a = 1.8 and f3 = 0.5, respectively. 

As the figure shows, data for the quantum regime is 
noisier than its classical counterpart. Finite size effects 
are also more prominent as can clearly be observed in the 
variation of the values of the slope 9. However, the trend 
observed is in agreement with the expected behavior in 
the thermodynamic limit. It suggests that the new quan- 
tum correlated critical points are indeed stable against 
the perturbation introduced by a short-range correlated 
potential. This opens up the prospect of measuring the 
new localization length exponents i^q in experiments in 
which random disorder with power-law correlations is in- 
troduced. The idea would be to engineer a random po- 
tential with desired properties which would then com- 
pete with the short-range disorder due to the presence 
of impurities in the semiconductor heterostructure. An 
example of such a system was described in Ref. where 
a thin magnetic film was placed in close proximity to 
the two-dimensional electron gas resulting in a random 
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FIG. 12: Quantum regime: slope 6 as a function of parameter 
It. K — indicates a purely long-range correlated potential, 
K = 1 corresponds to a purely short-range correlated poten- 
tial. 



magnetic field. By controlling the height fluctuations of 
the film one can in principle engineer a random potential 
with desired correlations. 



VII. CONCLUSIONS 

We have analyzed the effect of spatial correlations of 
the disorder potential on the localization-dclocalization 
transition in the integer quantum Hall system. The main 
purpose of the analysis was to understand separately the 
role played by quantum effects and disorder on the crit- 
ical properties of the transition by comparing classical 
and quantum localization. 

The analysis of electron localization in the classical 
case reveals that the localization length exponent is com- 
pletely determined by the statistical properties of the 
level lines of the random potential. Using the tools of 
percolation theory, we were able to implement well es- 
tablished numerical methods to perform a detailed check 
of the validity of the Harris criterion. We found numeri- 
cal support for the extended version of the criterion, con- 
firming previous theoretical arguments^. 

The quantum regime was analyzed by studying the real 
time density-density propagator as proposed in Refl2S 
This method, originally introduced as an alternative way 
to calculate the critical exponent v for a short-range 
correlated potential, proved to be an excellent testing 
ground to analyze the role of disorder. Wc obtained (fol- 
lowing a rather simplistic finite size analysis) a value for 
ly = 2.33 ± 0.09 remarkably close to experimentally mea- 
sured values and numerically calculated ones. We also 
found that when the decay of power-law correlations of 
the random potential is slow enough this can destabi- 
lize the quantum critical point and lead to new critical 
behavior, as predicted by the extended Harris criterion. 



Comparison between classical and quantum regimes 
with long-range correlated disorder potentials, indicate 
that the effect of disorder correlations are qualitatively 
similar for quantum and classical systems. The main dif- 
ference between these two cases seems to be the critical 
value of a (the exponent determining the long-range cor- 
relations) below which the value of i/ is changed. The 
value of a* ~ 0.85 is smaller compared to its classical 
counterpart a* = 1.5. This suggests that quantum fluc- 
tuations can be thought of as effectively smearing out the 
correlations in the random potential thus shifting a*. 

An immediate consequence of these results is that in 
the quantum case, there is no direct connection between 
the statistical properties of equipotential lines and the lo- 
calization length exponent. Specifically, we were able to 
show that the fractal geometry of the equipotential lines 
can be continuously varied while the quantum localiza- 
tion length exponent does not not change. 

The detailed numerical study of the variation of the 
critical exponent v with the power-law exponent a, led 
us to propose a precise statement for the quantum version 
of the extended Harris criterion. These numerical results 
are supported by a scaling argument, in close analogy 
with the classical case. 

In order to address the experimentally relevant situa- 
tion of samples with disorder potentials of different ori- 
gins and likely, different correlations, we analyzed the 
effect of competing correlations in disorder potentials on 
the IQHT. In the semi-classical limit, our numerical re- 
sults give support to previous theoretical arguments for 
the stability of the classical correlated-pcrcolation crit- 
ical points. In the quantum regime we observed that, 
while the value of the critical exponent v varies as the 
short-range correlated potential is introduced, the trend 
is similar to the one observed in the classical regime. 
As a consequence, our results show no qualitative differ- 
ence between classical and quantum regimes regarding 
the stability of these critical points. More importantly, 
they point to the possibility of seeing effects of the corre- 
lated quantum percolation critical points in experiments 
on two-dimensional electron gases in semiconductor het- 
erostructures where short-range disorder due to impuri- 
ties is unavoidable. 
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APPENDIX 

In this appendix we present a calculation based on per- 
turbation theory that shows that at short times the wave- 
packet has a ballistic spread, which is what was observed 
numerically. If \tp{t)) is the wave- function of the wave- 
packet at time t, the spread is given by: 
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imi^Am) = imi^^im) - mmmf (a.i) 

Here we analyze the first term in detail, while the sec- 
ond one can be dealt with in a similar way. 
In the Heisenberg representation: 

{m\^'\m) = {m\^'{t)\m) (a.2) 

where 

x'^{t)=e'"' x^{0)e-'"' . (A.3) 
Using the well known operator identity: 



imi^'im) imi^'mm) (a.4) 
+ it{m\[H;x'm\m) 

The probability distribution of energies is determined 
by the probability distribution of the disorder potential, 
which is chosen to be symmetric around the value V = 0. 
Since the wave-packet contains all the eigenstates, the 
averaging over disorder (average over energies) eliminates 
the linear term in the expansion above. Thus, the first 
non- vanishing contribution at short times is given by the 
quadratic term which corresponds to a ballistic spread. 
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